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Abstract 

We apply the coupled cluster method to high orders of approximation and exact diagonalizations 
to study the ground-state properties of the triangular-lattice spin-s Heisenberg antiferromagnet. 
We calculate the fundamental ground-state quantities, namely, the energy eo, the sublattice mag¬ 
netization Mgub) the in-plane spin stiffness ps and the in-plane magnetic susceptibility y for spin 
quantum numbers s = 1/2,1,..., Smax; where S m ax = 9/2 for eo and Smax = 4 for ps and 

•Smax = 3 for X- We use the data for s > 3/2 to estimate the leading quantum corrections to the 
classical values of eo, Mgub, Ps, and x- In addition, we study the magnetization process, the width 
of the 1/3 plateau as well as the sublattice magnetizations in the plateau state as a function of the 
spin quantum number s. 
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I. INTRODUCTION 


In the 1970s Anderson and Fazekas 


,Q 


first considered the quantum spin-1/2 Heisen¬ 
berg antiferromagnet (HAFM) for the geometrically frustrated triangular lattice and they 
proposed a liquid-like ground state (GS) without magnetic long-range order (LRO). Later 
on it was found that the spin-1/2 HAFM on the triangular lattice possesses semi-classical 


three-sublattice Neel order, see, e.g.. Refs. 


19| . However 


Msnh is drastically diminished in the s = 1/2 model 




he sublattice magnetization 


18| because of the interplay 


between quantum fluctuations and strong frustration. The small magnetic order parameter 


indicates that the semi-classical magnetic LRO is fragile and that smal 


the Hamiltonian may destroy the magnetic LRO, see, e.g.. Refs. 


20 


28|. 


additional terms in 


Although very precise data for the relevant GS quantities are available for unfrustrated 
HAFM’s on bipartite two-dimensional lattices, see, e.g.. Refs. 


29 


32| related to the square 


lattice, the corresponding data for the triangular lattice are less precise. This lack of precision 
is related to the strong frustration in the system that, e.g., does not allow one to apply 
the quantum Monte Garlo method. Moreover, the spin-wave approach is less efficient for 
frustrated lattices than it is for non-frustrated lattices. Nevertheless, spin-wave theories 
are considered as appropriate, in particular, if the spin quantum number s is not s = 
1/2 or s = 1. Perhaps the most accurate result for the GS order parameter (i.e., the 


sublattice magnetization Mg, 


for s = 1/2 has been obtained by a recent density matrix 


renormalization group study [16[, where a result of Msuh = 0.205 has been found. 


The continuous interest in the triangular-lattice HAFM is (last but not least) also re- 
ated to a fluctuation-induced magnetization plateau at 1/3 of the saturation magnetization 


33l-l53| . In particular, two model compounds, namely Ba 3 GoSb 209 with s = 1/2 and 


Ba 3 NiSb 209 with s = 1, have been shown very recently to demonstrate an excellent agree¬ 
ment between the experimentall y m easured magnetization curves and those curves from 


theoretical predictions, see Refs. (39|, l45|, l4^ for s = 1/2 and Refs. [4J, 


481 for s = 1. 


In the present paper we consider the Hamiltonian 


H = ( 1 ) 

(*i> * 

where the sum runs over nearest-neighbor bonds {ij) on the triangular lattice, (sj)^ = 
s(s -|- 1), and h is an external magnetic held. We consider arbitrary spin quantum number 


2 












































s. We use the coupled cluster method (CCM) to high orders of approximation to determine 
the GS properties in zero magnetic held, i.e., the GS energy per spin cq, the sublattice 
magnetization M^uh (order parameter), the spin stiffness ps, and the uniform susceptibility x- 
These quantities constitute the fundamental parameters determining the low-energy physics 
of the triangular Heisenberg antiferromagnet. Moreover, the stiffness and the susceptibility 
are used as input parameters in scaling functions for various observables 541. 

In addition to the zero-held quantities we also consider the magnetization process M{h) 
and determine the 1/3 plateau in the M(h)-curve. We complement the GGM calculations 
by carrying out Lanczos exact diagonalization of hnite lattices. 
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T 


^ III ^ 



FIG. 1: Reference states used for the CCM calculations. 


II. METHODS 


A. Lanczos exact diagonalization 


The Lanzcos exact diagonaliazion (ED) is one of the most useful methods that can be 


used to investigate frustrated quantum spin systems, see, e.g.. Refs, 
lattices of size iV = 36 are common for ED calculations for ^n s = 


^-62|. Although 


N accessible for ED shrinks signihcantly, see, e.g.. Refs. 


48 


57 


60 


63 


1% the system size 


-65|. Hence, we use 


the ED here in order to complement the results of the GGM (that yields results in the limit 
N —)■ oo). We use J. Schulenburg’s spmpacA: code[66| to calculate the magnetization curves 
for s = 1/2,1,..., 5/2. The maximum lattice size for s = 2 and 5/2 is N = 12, whereas for 
s = 3/2 we have results for N = 12,18, 21. For s = 1 the largest lattice we can consider is 
N = 27. We use these data to analyze the s-dependence of the 1/3 plateau. 
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B. Coupled cluster method 


The coupled cluster method (CCM) is a universal many-body method widely used in 
various fields of quantum many-body physics, see, e.g. Refs. j^, Meanwhile, the CCM 


has been established as an effective tool in the t 
see, e.g., the recent papers Q,y, 


39 


48 


69 


leory of frustrated quantum spin systems. 


80[ |. Here we illustrate only some features of 


the CCM relevant for the present paper. For more general information on the methodology 


of the CCM, see, e.g.. Refs. 


68 


81 


85|. 


The CCM calculation starts with the choice of a normalized reference state |<h). We choose 
the classical GS of the model as reference state, which is well known for the triangular HAFM 


for arbitrary fields, see, e.g.. Refs. 


35 


39 


40| and Fig. [T] For zero field it is three-sublattice 


Neel state, i.e., state I with a = 60° in Fig. [T] For finite magnetic fields non-collinear planar 
states with field dependent pitch angles a and {3 are classical GS’s, see Fig. [1] The reference 
state is a collinear state (so-called up-up-down state, see state II in Fig. [T]) only at the 1/3 
plateau. With respect to the corresponding reference state, we then define a set of mutually 
commuting multispin creation operators C/, which are themselves defined over a complete 
set of many-body configurations I. We perform a rotation of the local axis of the spins such 
that all spins in the reference state align along the negative z axis. The specific form of the 
spin-operator transformation depends on the pitch angles of the reference state. In this new 
set of local spin coordinates the reference state and the corresponding multispin creation 
operators C/ are given by 


|<l) = I = s+, s+s+ , s+s+4 ,... , 


( 2 ) 


where the indices n,m,k,... denote arbitrary lattice sites. In the rotated coordinate frame 
the Hamiltonian becomes dependent on the pitch angles. With the set {|$),C/} the CCM 
parametrization of the exact ket and bra GS eigenvectors |T) and (T| of the many-body 
system is given by 


|>I'>=£!®|4), S = 5^a,C+ 


/^O 


($1 = {4.|Se-y S = l + J]a,C,-, 

1^0 


( 3 ) 

( 4 ) 


where Cj = The CCM correlation operators, S and S, contain the correlation coef- 
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ficients, aj and a/, which can be determined by the CCM ket-state and bra-state equations 


($|C7e-^ife^|$) = 0 ; VJ ^ 0 
{<^\Se-^[H,C+]e^\<!>) = 0 ; VJ ^ 0. 


(5) 

( 6 ) 


Note that each ket-state equation belongs to a specific creation operator 
C/ = ■Sn'Sm) i-G., it correspouds to a specific set (configuration) of 

lattice sites n, m, fc,... . By using the Schrodinger equation, H\^) = we can 

write the GS energy as E = ($|e“'^ife‘^|$). The sublattice magnetization is given by 
Msnh = “(1/-^) 1^)’ where sf is expressed in the transformed coordinate system. 

The total magnetization M aligned in the direction of the applied magnetic field h in terms 
of the global axes prior to rotation of the local spin axes is given by M = {Ma + Mb + Mc)/3, 
where Mb, and Me are the magnetizations of the three individual sublattices, cf. Fig. [H 
given by 


M, 


A,B,C 


N 


A,B,C 




^A,B,C 


1 ^), 


(7) 


^A,B,C 


where the index runs over all Na sites on sublattice A, the index runs over all Nb 
sites on sublattice B, and the index ic runs over all Nc sites on sublattice C, and N = 
Na + Nb + Nc- The CCM results for the ground state energy and the total magnetization as 
a function of the magnetic field can be used to calculate the uniform magnetic susceptibility, 
given by 




dM 


1 d^E 


dh N dh^ 
Note that we consider here y as susceptibility per site 


( 8 ) 


The GS energy depends (in a certain CCM approximation, see below) on the pitch angles. 
In the quantum model the pitch angles may be different to the corresponding classical values. 
Therefore, we do not choose the classical result for the pitch angles in the quantum model. 
Indeed, we consider them as a free parameter in the CCM calculation, which has to be 
determined by minimization of the CCM GS energy with respect to the pitch angles. An 
exception is the zero-field case, where the pitch angle is fixed to a = 60° (the three-sublattice 
Neel state). 

The spin stiffness ps measures the increase of energy rotating the order parameter of 
a magnetically long-range ordered system along a given direction by a small twist (pitch) 


5 






angle 6 per unit length, i.e., 


E{9) E{e = 0) 1 


p,9^ + 0{9^), 


(9) 


N N 2 

where E{9) is the ground-state energy as a function of the twist angle. For the triangular 
lattice the twist is imposed along a attice basis vector and it is within the plane defined 
by the order parameter, see Refs, js, 10], where the twist along both directions leads to 


identical results jSTj j. 

For the many-body quantum system under consideration it is necessary to use approx¬ 
imation schemes in order to truncate the expansions of S and S in Eqs. ([3]) and (jl]) in a 

T^well established SUBn-n approximation scheme, cf., e.g., 
, where the correlation operators contain no more than n 


practical calculation. We use t 
Refs. Q, 27, 39, 48, 


70 


80 


82 


spin flips spanning a range of no more than n contiguous lattice sites ^ 


Using an efficient parallelized CCM code 89|| we are able to solve the CCM equations up 
to SUBlO-10 for s = 1/2 (where, e.g., for the zero-held case reference state a set of 1 054 841 
coupled ket-state equations has to be solved). For s > 1/2 the number of CCM equations 
increases noticeably. Hence, the highest order of approximation is then SUB8-8 (where, e.g., 
for the susceptibility for s = 3 a set of 2 179 007 equations has to be solved). 

The SUBn-n approximation becomes exact only for n oo. We extrapolate the ‘raw’ 
SUBn-n data to n —)■ oo. Much experience exists relating to the extrapolation of the 
GS energy per site eo(n) = E{n)/N, the magnetic order parameter Mguh, the spin stiff¬ 


ness ps, and the susceptibility y. Thus eo(n) = oq -l- ai + 


tested extrapolation ansatz for the GS energy per spin 


70 


73|, 


78 






/n)^ is a very well- 


85| . An appropriate 


extrapolation rule for the magnetic ord er p arameter of antiferomagnets with GS LRO is 
Msnh = bo + bi{l/n) + 52(l/n)2 3 , 1 ^ . [ 3 , 0 , 0 | • For the stiffness ps as well as for the 
susceptibility y we use the same rule as for M^uh, he., X{n) = cq + ci(l/n) -|- C 2 (l/n)^, 


X = Ps, X, which is a 
Ps well, see Refs. 


Die to describe the asymptotic behavior of the CCM-SUBn-n data for 
71], and y, see Ref. 39|. 


The selection of the SUBn-n data included in the extrapolation is a subtle issue. Often it 


is argued that the lowest-order data (i.e., SUB2-2 and SUB3-3) ought to be exc, 
the extrapolation because these points are rather far from the asymptotic regime 


uded from 


76 


78|. 


This argument is particularly valid for models which include larger-distance exchange bonds 
(e.g., so-called J 1 -J 2 models) j^, 71, 72 , 74 , 75|. However, for the triangular Heisenberg 
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antiferromagnet with only nearest-neighbor bonds even the lowest approximation orders £t 


well to the extrapolation 


Q,E, 


48| . Another point is the odd-even problem, i.e., for odd 


and even numbers n of the SUBn-n approximation the extrapolation may have different 


£t parameters 


73 


78| . However, this problem occurs primarily for bipartite systems (with 


collinear reference states, where no odd-numbered spin flips enter the correlation operators S 
and S'), whereas for noncollinear reference states (where odd-numbered spin flips are present 
in S and S) relevant for many frustrated systems both, odd and even SUBn-n, might be 


combined in one and the same extrapolation formula P,l76|,l78l. 

In order to £t the data to the extrapolation formulas given above (which contain three 
unknown parameters), it is desirable (as a rule) to have at least four data points to obtain a 
robust and stable £t. To obey this rule we apply here the above given ex trap olation formulas 
using (i) even n = 2, 4, 6, ... and (ii) odd and even n = 3, 4, 5, 6, 7, ... 90|. The maximum 
approximation level SUB?7,max-''^max for s = 1/2 is Umax = 10 for cq and M^uh and Umax = 9 
for y and whereas for s > 1/2 we have Umax = 8. In the following we call case (i) ’extral’ 
and case (ii) ’extra2’. The difference between both cases can be considered as a measure 
of accuracy of our CCM results. To illustrate our extrapolation procedure, we present in 
Fig. [2] the extrapolations of Cq, Msnh-, Xi and pa for s = 1/2 and s = 1. Obviously, the 
extrapolations work very well for eo, Msub, and ps, and, both schemes, ’extral’ and ’extra2’, 
lead to very similar results. There is some scattering of the SUBn-n data only for y, and, as 
a consequence, there is a visible difference between the extrapolations ’extral’ and ’extra2’. 


III. RESULTS 


The zero-field case 


In this section we present CCM results for the GS energy per spin eo, the sublattice 
magnetization (order parameter), the spin stiffness p^, and the uniform susceptibility 

y for spin quantum numbers s = 1/2,1,..., 9/2 (for Cq and Msub), for s = 1/2,1,..., 4 (for 
Pa ), and for s = 1/2,1,..., 3 (for y). Moreover, we use the data for s > 3/2 to estimate 
the leadin g q uantum corrections to the classical values to compare with the 1/s spin-wave 
expansion [Sl. la. la. Il7|. 

The data for the GS energy and the sublattice magnetization are collected in Table [H The 
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TABLE I: Extrapolated CCM results for the GS energy per spin, eoln^ooj the GS sublattice 
magnetization, M^uhln^oo, the spin stiffness, ps\n^oo, and the susceptibility, xln^oo- We mention 
that the spin-wave velocity Cgwt can be calculated from ps and y by using the hydrodynamic relation 
cLt = Ps/X- 





extral 



extra2 




eo/s^ 

Msub/S 

Psjs^ 

X 

eo/s^ 

Msub/S 

Ps/s^ 

X 

s 

= 1/2 

-2.2056 

0.4307 

0.3103 

0.0652 

-2.2045 

0.4248 

0.2990 

0.0553 

s 

= 1 

-1.8384 

0.7303 

0.6429 

0.0956 

-1.8367 

0.7350 

0.6572 

0.0902 

s 

= 3/2 

-1.7234 

0.8169 

0.7636 

0.0996 

-1.7223 

0.8232 

0.7757 

0.0972 

s 

= 2 

-1.6667 

0.8628 

0.8246 

0.1023 

-1.6659 

0.8695 

0.8342 

0.1012 

s 

= 5/2 

-1.6329 

0.8909 

0.8610 

0.1041 

-1.6323 

0.8973 

0.8687 

0.1035 

s 

= 3 

-1.6105 

0.9096 

0.8850 

0.1054 

-1.6000 

0.9155 

0.8913 

0.1050 

s 

= 7/2 

-1.5946 

0.9229 

0.9019 

- 

-1.5941 

0.9282 

0.9071 

- 

s 

= 4 

-1.5826 

0.9328 

0.9145 

- 

-1.5823 

0.9376 

0.9188 

- 

s 

= 9/2 

-1.5734 

0.9404 

- 

- 

-1.5731 

0.9448 

- 

- 


difference between both extrapolation schemes, ’extral’ and ’extra2’, is largest for lower spin 
quantum numbers s, although it is still small for all values of s. In the extreme quantum 
limit s = 1/2 the density matrix renormalization group result Q Msuh/s = 0.410 is slightly 
lower than our CCM result. 

Let us now compare our data for cq and Ms^b with recent higher-order spin-wave results 
by Chernyshev and Zhitomirsky 17|]. Chernyshev and Zhitomirsky found that eo(s) = 
-1.5s2(1 +0.218412/s + 0.0053525/s2) and Msub(s) = s(l-0.261303/s+ 0.0055225/5^). We 
£t our extrapolated CCM data for s = 3/2, 2,..., 9/2 using the ansatz 
X(s) = (1 — xi/s — X 2 /s^), X = eo, Msnh- The classical values are eo|s-).oo = —3s^/2, 

and MsubI s —>-00 

The values for the 1/s expansion parameters Xi and X 2 are listed in Table [ITl and the 
corresponding results are depicted in Fig. [3[ _ ^For the GS energy xi and X 2 are in very 


good agreement with the spin-wave results 


17j. The leading coefficient Xi for the order 














FIG. 2: Extrapolation of the CCM-SUBn-n for the GS energy eo (a), the sublattice magnetization 
-^^sub (b), the susceptibility y (c), and the spin stiffness ps (d) for spin quantum numbers s = 1/2 
and 1 using two different extrapolation schemes (labeled by ’extral’ and ’extra2’), cf. main text. 

parameter also fits well to the spin-wave term. However, we obtain (small) negative values 
instead of a (small) positive one for the next-order coefficient X 2 - The good agreement 
between the spin-wave and the CCM results for the GS energy is also evident in Fig. [3](a). 
Moreover, the 1/s expansion up to second order yields reasonable results for cq even for 
the extreme quantum case s = 1/2. On the other hand, the deviation for the sublattice 
magnetization becomes noticeable for s < 3/2, see Fig. [3](b). Thus, by contrast to the 
GS energy, the 1/s expansion of M^uh/s up to order s“^ leads to values for s = 1/2 with 
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FIG. 3; Extrapolated CCM data (schemes ’extral’ and ’extra2’, see main text) as a function 
of 1/s compared with higher-order spin-wave theory (SWT, blue line) taken from Ref.jl^. The 
symbols represent the CCM data points, the corresponding red and black lines show the fit function 
X(s) = X|s_^oo (l — x\ls — X 2 I where the data for s > 3/2 were used for the fit. (a): GS energy 
cq. (b) Sublattice magnetization Ms^b (the green symbol shows the density matrix renormalization 
group result (DMRG) for s = 1/2 from Ref. Q). 


limited accuracy. We know that the sublattice magnetization of the unfrustrated square- 


lattice s = 1/2 Heisenberg antiferromagne 


; obtained by higher-order spin-wave theory 301 


agrees well with quantum Monte Carlo 29|] and CCM|32|] results. Thus the deviation for 


the triangular lattice might be attributed to the enhanced quantum fluctuations caused by 
frustration leading to a particularly small value of for s = 1/2. 

Next we discuss the in-plane spin stiffness ps and the in-plane magnetic susceptibility y. 
Results are given in Table HI As already mentioned above, the extrapolation of the CCM- 
SUBn-n data works well for ps, although it is less accurate for x, i.e. for x the deviation 
between the schemes ’extral’ and ’extra2’ is noticeable, cf. Tabledand Fig. [2l The spin-wave 
large-s relations are ps = s^(l — 0.3392/s) (Ref. js]) and x = (1/9)(1 — 0.1425/s) (Ref. jol). 
Fitting the data from TableUfor s > 3/2 using the ansatz X{s) = (1 — Xi/s — X 2 / 

^ = Ps^Xi Ps|s->-oo = and x|s->-oo = 1/9, gives values for the 1/s expansion parameters 
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TABLE 11: Parameters xi and X 2 of the 1/s expansion X(s) = (l — xi/s — X 2 /s^) obtained 

from the extrapolated CCM results for the GS energy eg, the sublattice magnetization Msuh, the 
spin stiffness ps and the susceptibility y. 



eo 

T/sub 



Xl 

X2 

Xl 

X2 

extral 

0.2176 

0.0071 

-0.2671 

-0.0120 

extra2 

0.2186 

0.0073 

-0.2416 

-0.0362 


Ps 

X 



Xl 

X2 

Xl 

X2 

extral 

-0.3355 

-0.0292 

-0.1587 

0.0045 

extra2 

-0.3166 

-0.0297 

-0.1440 

-0.0662 


Xi and X 2 which are listed in Table [TTl The leading coefficient xi hts reasonably well to the 
spin-wave term. The deviation between the schemes ’extral’ and ’extra2’ for y results in 
different signs of the second-order term X 2 . We show the 1/s dependence of ps and y in 
Fig. m It is obvious that the large-s approach works surprisingly well for ps down to s = 1/2, 
whereas it seems to fail for x for s = 1/2. 


B. The magnetization process 


As already mentioned in the introduction, t 
was investigated previously in numerous papers 


le magne 


34-36 


38 


iza 


39 


ion process M(h) for s = 1/2 


M{h) curve for the ^ecihc case s = 1 was much less studied 


4 


42 


45 


m 


49 


51 


The 


large-s approaches 


35 


40 


42 


50 


60| . Several quasi-classical 


52| can be used to obtain an estimate for M{h) and the 


1/3 plateau also for lower spin quantum numbers. However, it is likely that these results 
have limited accuracy (see also the discussion in Sec. IIII A|) . We mention again that very 
good agreement between experimental and theoretical CCM data for s = 1/2 and s = 1 has 


been reported 


45 


48| very recently. 


Let us mention that the CCM calculations of the M{h) curves are extremely time con- 
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FIG. 4: Extrapolated COM data (schemes ’extral’ and ’extra2’, see main text) for the spin stiffness 
Ps (a) and the susceptibility x (b) as a function of 1/s compared with spin-wave theory (SWT, blue 
line) taken from Ref. SJ (for ps) and Ref. ^ (for x)- The symbols represent the COM data points, 
the corresponding red and black lines show the fit function 2f(s) = X|s_).oo (l — xi/s — X 2 /s^), 
where for the fit the data for s > 3/2 were used. 


suming, because the field dependent quantum pitch angles for each value of the magnetic 
held have to be determined by minimization of the CCM-SUBn-n GS energy with respect 
to the pitch angles. Hence we consider only even SUBn-n approximations until n = 6. We 
have data for SUB8-8 only for the critical helds, hd and hc 2 , which bound the 1/3 plateau, 
and only for the most interesting extreme quantum cases of s = 1 and s = 1/2. We show 
the CCM-SUB6-6 magnetization curves for s = 1/2,1,3/2, 2, and 5/2 in the main panel of 
Fig. [5l The width of the 1/3 plateau shrinks with increasing spin quantum number s which 
is clearly seen in Fig. [5l From the experimental point of view this shrinking of the plateau 
width is relevant. Thus for s = 2 the plateau width {hc 2 — hd)/s is only about 25% of the 
width for s = 1/2 which makes its detection for large s by measurements at low (but hnite) 
temperatures more challenging. From Fig. [5] it is obvious that all curves for s > 1/2 are 
close to each other. Below and above the plateau they show almost the classical linear h 
dependence of M. The s = 1/2 curve is well separated and shows a pronounced deviation 
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FIG. 5: Main: The magnetization curve calculated within CCM-SUB6-6 approximation for 
s = 1/2,1,3/2, 2 and 5/2. Inset: The magnetization curve for s = 3/2 calculated with exact 
diagonalization for finite lattices of = 18 and 21 sites compared with the CCM-SUB6-6 data. 

from linearity. The curves shown in the inset demonstrate that CCM and ED data agree 

well. 

The CCM approach allows to calculate the individual sublattice magnetizations Ma = 
Mb and Me as well as the quantum pitch angles a and /9, cf. Fig. [T] These quantities, are 
accessible, in principle, in neutron scattering experiments. They provide a deeper insight 
in the details of the magnetization process and the role of quantum fluctuations. We show 
Ma, Mb and Me in Fig. |6]and a and (3 in Fig. [71 An interesting feature is the non-mo notonic 
behavior of Me for h < hd (present only in the quantum model) and of Ma, Mb and a 
above the plateau. There is a strong increase in the slopes (except for Me and (3 near hd) 
as one one approaches the plateau from below or above. 

For the collinear plateau state at one third of the saturation (the so-called ’up-up-down’ 
state, see state II in Fig.[T]) we have calculated SUBn-n data of the sublattice magnetizations 
Mup = Ma = Mb and Mdown = Me up to rimax = 10 for s = 1/2 and 1 and up to rimax = 8 
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FIG. 6: Field dependence of the sublattice magnetizations = Mb (left) and Me (right), cf. 
Fig. [H calculated within CCM-SUB6-6 approximation for s = 1/2,1, 3/2, 2 and 5/2. 



FIG. 7: Field dependence of the pitch angles a (left) and /3 (right), cf. Fig. [H calculated within 
CCM-SUB6-6 approximation for s = 1/2,1,3/2, 2 and 5/2. 

for s = 3/2, 2,..., 4. Again we perform an extrapolation of the CCM-SUBn-n data applying 
Mx = bo + bi{l/n) + b 2 {l/n)‘^, X = A,B,C. We use only even SUBn-n approximations 
for the extrapolation and the corresponding scheme is therefore ’extral’ (see Sec. Ill Bl) . The 
resulting data for = Mb and Mq are given in Table HTTl We hnd that Mup is always larger 
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than |Mdown|- In particular, the difference in the magnitude of both quantities is remarkably 
large for s = 1/2. As for the zero-field case, the sublattice magnetizations within the plateau 
state are reduced by quantum fluctuations. This reduction is, however, much smaller than 
that for the canted zero-field state. 

We obtain the 1/s dependence for Ma = Mb and Me by htting our extrapolated CCM 
data for s > 3/2 using (as previously) the ansatz vA(s) = X|s_).oo (1 — xi/s — X 2 /s^), X = 
Ma,b,c- The classical values are Ma,b\s^oo = s and Mc\s^oo = —s- The values for 1/s 
expansion parameters xi and X2 are xi = —0.1003 and X2 = 0.0091 for Ma,b and xi = 
-1-0.2006 and X 2 = —0.0182 for Me and the corresponding 1/s behavior is shown in Fig. [HI 

TABLE III: Extrapolated CCM results (extrapolation scheme ’extral’) for the sublattice magne¬ 
tizations in the plateau state, Mx\n^oo-, X = A,B,C. 


s = 1/2 

s = l s = 3/2 

s = 2 s = 5/2 

s = 3 s = 7/2 

s = 4 

Ma,b 0.8392 

0.9095 

0.9372 

0.9521 

0.9614 

0.9676 

0.9721 

0.9755 

Me -0.6783 

-0.8190 

-0.8743 

-0.9043 

-0.9227 

-0.9352 

-0.9441 

-0.9509 


We shall now discuss in some detail the critical fields, hd and hc 2 , that define the position 
and the width of the plateau in the M{h) curve. These results are also relevant for the 


experimental searches of the mage 


ization plateau in magnetic compounds with triangular- 


lattice structure, cf. Refs. j44j-l46l. l48[| . We present our data in Table HVl and show the 1/s 
dependences of and hc 2 in Fig. [91 In addition to our CCM and ED data, we show also 
relevant data for hd and hc 2 from Refs, j^, Q, [^, for comparison. The monotonic 
shrinking of the plateau widths known for the large-s approaches is also clearly seen in our 
CCM and ED data. 

We notice that the plateau ends at hd and hc 2 behave differently with increasing s. 
Although the lower plateau end is only slightly shifted, the shift of the upper one at hc 2 is 
more pronounced, see also Figs.[5}[Zl We also see that the data for hd and hc 2 provided in the 
literature for the extreme quantum case s = 1/2 exhibit a rather large amount of scattering. 
As already found for the zero-field sublattice magnetization and susceptibility, cf. Sec. IIII Al 


the quasi-classical large-s approaches |35|, [^, [52| for s = 1/2 noticeably deviate from our 


data directly calculated for s = 1/2. Moreover, the recent real-space perturbation theory|52|] 
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FIG. 8: Extrapolated CCM data (scheme ’extral’, see main text) for the sublattice mag¬ 

netizations in the plateau state, Ma = Mb and Mq, as a function of 1/s. The symbols 
represent the CCM data points, the corresponding red and black lines show the fit function 
Mx{s) = Mxls^oo (1 -xi/s- X 2 /s^). 


yields an 1/s-dependence of hd and hc 2 that significantly deviates from the spin-wave, ED 
as well as the CCM behavior. On the other hand, the recent large-scale cluster mean-field 


approach of Ref. 


5l[ | yields hd/s = 2.690 and h^/s = 4.226 for the s = 1/2 case, and these 


values are close to the CCM-SUB8-8 results. 


IV. SUMMARY 


The HAEM on the triangular lattice is a basic model of quantum magnetism. The 
theoretical treatment of this frustrated spin model is challenging and the precision of the 
existing data for the basic parameters determining the low-energy physics of the model 


is less than that of corresponding unfrustrated models such as the square- 


Furthermore several magnetic compounds exist, see e.g. Refs. 


44 


46 


91 


attice HAEM. 


94| . that are 


described by this model. Hence, there is a need to improve the accuracy of the data available 
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TABLE IV: 
s = 1/2 ED 


Critical fields 
data also for N = 


and hc 2 , where the one-third plateau begins and ends. Note that for 
: 36 and 39 are available 3g, Ssj, 61]. 


CCM 


SUB2-2 

SUB4-4 

SUB6-6 

SUB8-8 

s 

hells 

he2ls 

hells 

hc2ls 

hells 

hc2/s 

hells 

he2ls 

1/2 

2.382 

4.344 

2.624 

4.482 

2.714 

4.370 

2.740 

4.290 

1 

2.727 

3.617 

2.788 

3.674 

2.809 

3.648 

2.814 

3.637 

3/2 

2.826 

3.405 

2.854 

3.440 

2.862 

3.429 

- 

- 

2 

2.873 

3.303 

2.888 

3.327 

2.892 

3.322 

- 

- 

5/2 

2.900 

3.242 

2.909 

3.261 

2.911 

3.257 

- 

- 

3 

2.917 

3.202 

2.923 

3.217 

2.925 

3.214 

- 

- 

ED 


N = 

12 

N = 

18 

N = 

21 

A = 

27 

s 

hells 

he2ls 

hells 

hc2ls 

hells 

hc2/s 

hells 

he2ls 

1/2 

2.615 

4.934 

2.805 

4.578 

2.794 

4.425 

2.745 

4.382 

1 

2.816 

3.916 

2.879 

3.786 

2.851 

3.717 

2.817 

3.695 

3/2 

2.864 

3.613 

2.913 

3.523 

2.890 

3.479 

- 

- 

2 

2.893 

3.461 

2.932 

3.391 

- 

- 

- 

- 

5/2 

2.911 

3.369 

- 


- 


- 

- 

3 

2.924 

3.308 

- 

- 

- 

- 

- 

- 


from theoretical investigations also from the experimental point of view. 

In the present paper we present large-scale CCM calculations for the basic GS parameters, 
energy cq, sublattice magnetization in-plane spin stiffness ps and in-plane magnetic 

susceptibility x, for arbitrary spin quantum number s. In addition to these zero-held quan¬ 
tities, we also consider the magnetization process. It is known from many previous studies 
for other frustrated quantum spin models, such as the kagome HAFM and the J 1 -J 2 square- 
lattice HAFM, that the CCM provides accurate GS results. Hence, the results presented 
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FIG. 9: Critical fields, hd and hc 2 , that bound the 1/3 plateau versus the inverse of the spin 


quantum number, 1/s, calculated by exact diagonalization for = 21 and 27 (labe’ 


ED 27) and CCM (labels SUB4-4 and SUB6-6) compared to large-s approaches 3a, l42 


cluster mean-field data 


s ED 21 and 


and 


5ll |. The blue lines (labeled by MZ) correspond to hd/s = 3 — l/6s^ and 


hc 2 /s = 3 -|- l/3s -|- l/6s^ (Ref. 5^) and the red lines (labeled by CG) correspond to hd/s = 


3 — 0.028/s and hc 2 /s = 3 


from Ref. and Ref. 


0.0717/s (Ref. 35]). The abbreviations YMD and TTZ represent data 


42l |. respectively. 


here will contribute to improve available theoretical data, especially for the extreme quan¬ 
tum cases s = 1/2 and s = 1 where large-s (spin-wave) theories do not to provide sufficently 
precise data. 

Although the present study is purely theoretical, the data presented here might be used to 
interpret experimental results for ma gnet ic compounds that are described by the triangular- 


lattice HAFM, see e.g. Refs. [44 


46 


91 


-94j. As already mentioned above, the compounds 
Ba 3 CoSb 209 and Ba 3 NiSb 209 are described well by the Heisenberg model considered here 


with spin quantum number s = 1/2 (Ba 3 CoSb 209 ) 45] and s = 1 (Ba 3 NiSb 209 ) 4J,l48|. We 
remark that very good agreement between the theoretical CCM data and experimental data 
has been reported for these cases. The CCM data presented here for s = 3/2 and s = 2 
might be useful for further studies of the magnetic compounds La 2 Ca 2 Mn 07 (with s = 3/2) 
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93|, l94| and FeGa 2 S 4 (with s = 2) 9ll. l92|. 
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